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We show that nonradiative interactions between atomic dipoles placed in a waveguide can give 
rise to deterministic entanglement at ranges much larger than their resonant wavelength. The 
range increases as the dipole-resonance approaches the waveguide's cutoff frequency, caused by the 
giant density of photon modes near cutoff, a regime where the standard (perturbative) Markov 
approximation fails. We provide analytical theories for both the Markovian and non-Markovian 
regimes, supported by numerical simulations, and discuss possible experimental realizations. 
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Introduction. — Dipoles can interact via photon ex- 
change, resulting in excitation transfer and mutual en- 
tanglement pp. When the interaction is mediated by 
radiation, i.e. real photons, it constitutes a dissipa- 
tive and hence quantum-mechanically incoherent pro- 
cess, whereby the generation of entanglement is gener- 
ally probabilistic [21 [3] , although certain entangled states 
are deterministically obtainable by engineering/control 
of the bath [3]. In this study, we are concerned with the 
nonradiative interaction that stems from the collective 
coupling of atomic dipoles to a common "bath" of pho- 
tonic modes [S]. Such nonradiative (dispersive) interac- 
tions are possible via their near or evanescent fields [5J. 
Quantum mechanically they are described as exchange 
of virtual, i.e. non-resonant, photons between the atoms, 
known as resonant dipole-dipole interaction (RDDI) 0- 
9 . In free space RDDI is dominant over radiation only at 
distances shorter than the resonant wavelength. Here we 
predict modified RDDI along with suppressed radiation 
in confined geometries, giving rise to coherent interaction 
at distances much longer than the resonant wavelength. 
This constitutes a novel route towards high-fidelity long- 
range deterministic entanglement. The principle that al- 
lows to appropriately modify the radiative and dispersive 
interactions is that they are mediated by the geometry- 
dependent field modes, populated by either real or virtual 
photons, respectively. Hence, the distance-dependence 
of the interactions is determined by the geometry. For 
example, when mediated by surface-plasmon-polariton 
modes in one dimension, both interactions appear to have 
long-range character, yet they are hindered by dissipation 
mechanisms [T0TU2| . E.g., in [TU], the radiative interac- 
tion sets the bound of the concurrence (entanglement) at 
C = 0.5. This bound is circumvented by a promising ap- 
proach to a coherent phase-gate based on the difference 
between super- and sub-radiant decay rates |llj . Still, 
ohmic losses and radiation to free-space modes may prac- 
tically limit the phase-gate operation to distances smaller 
than a wavelength. Radiation, however, can be sup- 
pressed in geometries that create cutoffs or bandgaps in 
the photonic mode spectrum. In such geometries RDDI 
can be drastically modified [TMT5"] . 

In our approach, photonic cutoffs or bandgaps are used 
not only to suppress radiation but also to enhance RDDI 



so as to make it the dominant effect. Our main result, 
obtained by essentially exact (nonperturbative) calcula- 
tions, is the possibility of extremely long-distance RDDI 
almost without radiation, and correspondingly high con- 
currence (nearly-perfect entanglement). This effect is 
predicted in waveguides for pairs of atoms whose dipo- 
lar transition frequency is just below the cutoff or band- 
edge of the waveguide. We thereby reveal the key princi- 
ple that enables coherent long-range interaction, poten- 
tially much stronger than possible decoherence effects, 
namely, the very large density of photon states near the 
cutoff. Thus, the enhancement of density of states due 
to the cutoff is reminiscent of that obtained using a cav- 
ity. However, unlike a cavity, the waveguide geometry 
is open along the propagation axis and does not restrict 
the separation of the atoms. In the Markov approxima- 
tion, the RDDI diminishes with the interatomic distance 
z as e~~ z /£, where £ increases as the atomic frequency 
approaches the cutoff (bandedge), allowing for entangle- 
ment at long distances. Yet, the standard Markov ap- 
proximation fails close to cutoff, which requires a nonper- 
turbative analysis, supported by numerical calculations. 

The model. — We consider a pair of atoms, modeled by 
identical two-level-systems (TLS) with energy levels \g) 
and |e) and transition frequency uj a . These are coupled to 
the vacuum field of a non-leaky waveguide, i.e. we neglect 
the TLS coupling to modes outside the waveguide - a rel- 
evant assumption in the situation considered below. The 

TLS-ficld dipole couplings are gta — \J ' u fc( r a), *a 
being the location of atom a = 1, 2, d the dipole matrix 
element of the \g) <-> \e) transition (taken to be real), 
and ujk and Ufc(r) the fc'th mode frequency and spatial 
function. The corresponding Hamiltonian in the dipole 
approximation 116) reads, in the interaction picture, 
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H AF = h^Y, [ l »«^ e "" fct + h - c -] [^" fa,< + h - c -] > 

a— 1 k 

(i) 

afc, a a being the mode and the TLS lowering opera- 
tors, respectively. In what follows, we analyze the atomic 
dynamics under the perturbative Markov approximation 
and without it. 

Markovian theory. — Adopting an open-system ap- 
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FIG. 1: (color online), (a) Schematic picture of our model, 
(b) Dynamics of excited state population for the atoms af- 
fected by the TMn mode well below cutoff, with [Eqs. (6)- 
(9)] wn = 500rn, uj a - u>u = -lOOrn and z 12 = 0.5\ a : 
Markovian theory (t) and simulation (s) results, (c) Diver- 
gence near cutoff uu of the single-atom spectrum G aa for 
TMn mode, Eq. jjj). 



7aa', leaving intact the coherent dynamics governed by 
Hdd in Eq. (|3|. Then, if initially only atom 1 is excited, 
we get a periodic exchange of the excitation between the 
atoms, at a rate A12, in the two-atom state 

\ipi2(t)) =cos(Ai 2 i)|ei,flb) +isin(Ai 2 t)|5i,e 2 ), (5) 

that superposes singly-excited product states of atoms 1 
and 2. A maximally entangled state is then achieved at 
odd multiples of the time t = tt/(4A 12 ). 

In order to illustrate how the radiative effects 7 QCt < 
can be suppressed we first consider the case of atoms 
placed inside a rectangular hollow metallic waveguide 
(MWG), with longitudinal axis z and transverse dimen- 
sions a and b (see Fig. la). Nonideal MWG and opti- 
cal fiber realizations will be addressed below. The atom 
interacts only with the MWG field modes TE mn (trans- 
verse electric) and TM mn (transverse magnetic) labeled 
by non-negative integers m,n |20j (see Appendix). Each 
TE/TM mn transverse mode has its own cutoff frequency 
Lu mn and dispersion relation w™", k z being the longitu- 
dinal wavenumber, 
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proach for the problem [17) , we identify the two atoms 
as the system and the continuum of EM vacuum modes 
as a bath, and consider the effects of the bath on the 
system (Fig. 1(a)). These are dissipative and dispersive 
effects that are related by the Kramers-Kronig relation 
and determined by the bath's two-point (autocorrelation) 
spectrum G aa '(ui), defined via 



^ 9k a g* ka > 

k 



dloG aa ' (lj). 



(2) 



From Fermi's Golden rule we obtain the rate of dissi- 
pation by radiation j aa i = 2irG aa ' (co a ), which for a = a 1 
represents the single-atom spontaneous emission rate to 
the guided modes and for a ^ a' describes the two-atom, 
distance-dependent, cooperative emission [18]. The dis- 
persive effect is obtained by second-order perturbation 
theory for the energy correction (cooperative Lamb shift 
[19]) of the two-atom states, associated with the bath- 
induced dipole-dipole Hamiltonian term, 



Hdd = -h- 4 o«' {vi& a > + & a *a>) , 



where A r 



Aa a ) r + A a a',+ and 



A QQ ' j= p — P 



du) 



G aa ' (oj) 



(3) 



(4) 



P denoting the principal value. The dissipative, inco- 
herent effect of 7 QQ ' gives rise to probabilistic interac- 
tion between the atoms. Hence, in order to achieve non- 
radiative, deterministic interaction we need a vanishing 



where uj^ — wjL"™ is the frequency of the k — 
TE /TM mnt k z mode, and c is the speed of light. The 
contribution of a specific transverse mode X mn (A = 
TE,TM) to the bath spectrum in Eq. S is obtained 
from the dispersion relation k z (ui) (Eq. (161) upon iden- 
tifying uj™ 1 — u, 
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Q(x) being the Heaviside step function. At this stage two 
key features of the waveguide structure must be noted: 
(1) below the cutoff io mn no \ rnn guided photon modes 
exist, and (2) the density of states ^* diverges near the 
cutoff. In what follows, we use feature (1) to suppress ra- 
diation and feature (2) to obtain long-distance and strong 
RDDI. 

In order to facilitate the analysis it is sufficient to con- 
sider the case where the atoms are polarizable only in 
the z direction, d = d z e z (for other polarizations see the 
Appendix or Ref. |15j). Since TE modes have a van- 
ishing z component of the electric field, only TM modes 
contribute to the bath spectrum, 
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i — i is introduced, where 

7T€Qacab 

d z sin(^a; a )sin(^y Q ) and x a ,y a is the 
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transverse position of atom a. Also note that k z is a 
function of w by virtue of Eq. ^ . 

Now, consider the case where the atomic resonance 
is below the lowest cutoff frequency, w a < w n for TM 
modes. Then, the atomic dipoles are not resonant with 
any of the field modes and radiation is suppressed, 7 QQ ' — 
2nG aa > (uj a ) = 0, from Eq. ((9 ). We are thus left only with 
the nonradiative RDDI Eq. |4b, 



A12 = 

where z\ 2 = 
is 



r, 



t 



y/i - (io a /uj mn )' 



-_e (r, 



(10) 



\zi — z 2 \ and the effective interaction range 



u mn a/1 - (u> a /uj mn y 



(11) 



These Markovian-theory results [H] predict that radia- 
tive dissipation is absent, while the RDDI decays expo- 
nentially with interatomic distance, typical of interaction 
mediated by evanescent waves. Yet, remarkably, Eqs. 
( 10 1 and (111 imply that as the atomic resonance w ap- 



proaches the lowest cutoff uin from below, the RDDI di- 
verges owing to the contribution of the TM\\ mode, and 
so does its range, determined by £n. This potentially 
enables deterministic generation of entanglement at very 
large distances. 

In order to test the above results we performed direct 
numerical simulations of the Schrodinger equation for the 
Hamiltonian ([I]), taking only the dominant m = 1, n = 1 
mode into account (see Appendix). Fig. 1(b) portraits 
typical dynamics of the atoms' populations, along with 
their entanglement, quantified by the concurrence |21j . 
for z±2 = 0.5A a with A a the atomic transition wavelength. 
As expected from Eq. ([5| , the maximally entangled state 
is generated at half the oscillation period of the popula- 
tion exchange. It is also apparent that when uj a is not 
too close to the cutoff uin, the simulation results agree 
with those of the Markovian analysis, Eqs. ^ and (10), 
within numerical accuracy. 

Validity of the Markovian theory. — The Markov ap- 
proximation used above breaks down as uj a approaches 
the cutoff. The general conditions for the validity of the 
Markov approximation reduce in our case to (see Ap- 
pendix) 



A 12 (u) a )A'l 2 (u a ) < 1, 



(12) 



where A 12 (w a ) is given by Eq. ( 10 ) and A'{ 2 (u} a ) is its sec- 
ond derivative w.r.t w a . In the limit u a — > ui mn , Ai2(w a ) 
and A' 1 ' 2 (wa) become singular and condition \12\ is not 
satisfied, as seen from Eq. (10) and Fig. 1(c). Thus, a 



non-Markovian theory is required in order to fully ana- 
lyze the possibility of long-distance RDDI and entangle- 
ment. Non-Markovian analysis has been performed be- 
fore for a single atom coupled to a continuum with a cut- 
off [551 H3j ; yielding the possibility of incomplete decay: 
decay of the excited state population to a steady-state 



value different from zero, as a result of the formation of 
atom-photon bound-states. Nevertheless, the Markovian 
analysis is very useful for RDDI in cases where nearly- 
complete entanglement (e.g. C > 0.95) is to be achieved, 
as seen below. 

Non-Markovian theory. — In order to account for the 
situation where uj a approaches the cutoff, we develop a 
nonperturbative and non-Markovian theory for RDDI, in 
the spirit of jSHj. From Hamiltonian ([I]), assuming that 
only atom 1 is initially excited, the state of the com- 
bined (atoms+modes) system can be written within the 
rotating-wave-approximation |16j as 

|^(*)) = oi (t)\e x , 5 2, 0) +a 2 (t) \ 9l , e 2 , 0) +^ b k (t) \ 9l , g 2 , l k ) 

k 

Inserting this state into the Schrodinger equation, we 
obtain dynamical equations for ai(t), a 2 {t) and bk(t). 
As before, we consider only the MWG transverse mode 
m = l,n = 1, this time for uj a close to the cutoff u>n, 
such that the denominator of the spect rum Q is w ell ap- 
proximated by \J (w/wn) 2 — 1 ~ \/2y^uj/ujii — 1. Using 
the Laplace transform in order to solve the dynamical 
equations, we then obtain the dynamics of the first atom 
(more details can be found in the Appendix), 



ai{t) = Vie- lUllt ^2< 
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Here 
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the roots 



k(-ViUjVt) 
(13) 

of d(u) = u 5 + 2W a u 3 - 
Wn, Cj = n(uj)/d'(uj) with 



^75 J- nv^ii" 
where W a = 

n(u) = -i(u 3 + W a u - 2^r Uv /5n), and F(u) 
(g-2 Z i 2 (Vwir/c)« _ where F(u) is expanded in Taylor 
series up to 5th order in u (Appendix) . The conditions of 
validity for this theory are thus given by the approxima- 
tion of the spectrum and the expansion of F(u), yielding 
< i an d zi 2 < (D'^f^./ ,, ^ , , respec- 

4un ^ V 4 1 2-rr u\\ Ay 2 \u\\— uj a \ ' " 

tively. However, in practice, another limitation on the 
precision of the theory comes from the numerical calcu- 
lation of the roots of d(u). 

Fig. 2(a) presents the dynamics of the atomic 
populations and interatomic entanglement in the non- 
Markovian regime. Very good agreement between the 
above theory and numerical simulations is observed. The 
main feature of the dynamics are Rabi-like oscillations 
similar to those of the Markovian case. Nevertheless, 
their amplitude is decreased as a result of excitation loss 
to the field modes by incomplete decay, setting the up- 
per bound on the achievable entanglement. Hence, as u> a 
approaches the cutoff, while the inter-atomic distance z 
is kept fixed, we get a tradeoff between increased RDDI 
strength and decreased maximum entanglement. This is 
shown in Fig. 2(b), where uj a is varied from very far 
from cutoff, where Markovian theory predictions Ai 2j m 
and C max = 1 apply, to very close to the cutoff, where 
A12 increases on the expense of C max . 
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FIG. 2: (color online). Atomic dynamics affected by TM\\ 
mode in the non-Markovian regime, (a) Atomic excitation 
probabilities and concurrence as a function of time. Here 
ojii = 500rn, 2i2 = A a and u a — un = — lOrn. The simula- 
tions (s) well agree with the theory (t), Eq. ( 13 1, (b) Tradeoff 



between RDDI strength A12 and maximal achievable concur- 
rence Cmax as a function of the atomic-resonance mismatch 
with cutoff (simulation results for uu = 500Fii, 212 = 0.5A a ), 
compared with Markovian theory Ai2,m, Eq. (10 1. The esti- 



mates for A12 are extracted by fitting the dynamics of simula- 
tions for various uj a values to Eq. |5| . (c) Entanglement gen- 
eration as a function of time at long-distance: the excitation of 
the first atom is calculated by Eq. ( |13[ ) and the concurrence is 
bounded by the maximal value of the plot. MWG realization 
(blue thin line): u>u = 2.17x 10 10 rn, u a -u n = -2xl0 4 Fn, 
z = 100 X a . Fiber-Bragg-grating realization (red thick line): 
wu = 6 x 10 7 rn, w»-wu = -1500rn, z = 20A a . (d) Fiber- 
Bragg-grating scheme: the atoms are coupled to the guided 
modes by their transverse evanescent tails [28] or when in- 
serted into the fiber 1271. 



Long-distance entanglement and possible 
realizations. — Using the analytical theory above, 
we shall now illustrate the possibility of long-distance 
entanglement by two examples. First, consider Rydberg 
atoms that pass through a cold metallic waveguide 
(MWG), similar to the setup in [Ml Hi] where the 
MWG replaces the superconducting cavity. The states 
\g) and |e) are the two circular states with principal 
quantum numbers 51 and 50, with transition frequency 
and dipole moment uj a — 2ir x 51.1GHz and d ~ 1250eao 
respectively, with e the charge of an electron and ao 
the Bohr radius [25]. Near the cutoff, Tn is similar 
to the free-space |e) — > \g) decay rate, estimated to 



z = 100A Q are plotted in Fig. 2(c), where A a ~ 6mm 
is the atomic wavelength, such that we obtain entan- 
glement with concurrence C = 0.983, at a distance 
z ~ 0.6m and for interaction time t 0.2ms [Fig. 
2(c)]. Considering possible imperfections we derived 
the dissipation rate due to ohmic losses of the atom- 
induced evanescent fields in a square waveguide (a = 6), 
"/loss < l^ii w ith Mo the vacuum permeability and R s 
the surface resistance (see Appendix). Normal metals 
may limit the achievable entanglement distance and 
fidelity as in [llj . However, for niobium superconducting 
plates at temperature T < IK, we take, as in [24], 
R s = 75nf2, yielding, for a w 6mm, 7; oss = 19.89Hz, 
much slower than the 0.2ms required for entanglement. 
Such a temperature also ensures that the thermal photon 
occupancy at u> a is negligible. In addition, as analyzed 
in [5S] , surface roughness of the metal plate may slightly 
change the mode structure and the location of the cutoff 
frequency, and correspondingly the calculated RDDI 
rate. Nevertheless, a cutoff below which the modes 
become evanescent with diverging density of states 
persists, hence the principle of our scheme still applies. 
Regarding our initial assumption of isolated waveguide 
modes, we recall that uj a is much smaller than the 
typical plasma frequency in metals (~ 10 16 Hz), so that 
the isolated modes of a perfect-conductor used here, are 
indeed adequate. 

Another option is that of optical fiber modes coupled 
to the atoms [57J[55]. Although the fiber's guided modes 
also possess cutoffs, they lack the two important features 
that we have highlighted for the MWG: (1) below cutoff 
the atoms are coupled to outside modes, hence sponta- 
neous emission exists at a rate comparable to that in 
free space; (2) the group velocity does not vanish 



at the fiber cutoff so that the density of states j^j- does 
not diverge. We can restore the second feature by con- 
sidering a fiber-Bragg-grating [53] : then, for a transverse 
fiber mode with dispersion u)(k z ), the group velocity does 
vanish at the bandedge of the uj spectrum corresponding 
to k z — ir /(An), with A the period of the grating and 
n the average refractive index [Fig. 2(d)]. The disper- 
sion near the upper boundary of the gap, u u , can be ap- 
proximately written as uj ss lo u + B(k z — ir/(An)) 2 with 
constant B, so that oc l/y'w — lj u diverges at uj u 
in the same way assumed in our non-Markovian theory 
(see Appendix for more details). Then, the atom can 
still emit to outside modes, but just below the bandedge 
oj u , RDDI, which is mediated by evanescent waves in the 
gap, can become much stronger and more long-distance, 
due to the divergence. We consider optical atomic tran- 
sitions, e.g. the D2 line of 87 Rb atoms with A a ~ 780nm 
and natural linewidth 2-7r x 6.07MHz. The results for 
z = 20A a ~ 16/zm are plotted in Fig. 2(c), yielding con- 
currence C = 0.9605 after t s» 3.55ns of interaction. 

Conclusions. — To conclude, the main result of this 
study is the demonstration of the possibility of long- 
distance interaction between dipoles by a nonradiative, 
deterministic and coherent process (RDDI) that is cru- 
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14.7Hz. The corresponding dynamics for 
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cially dependent on the waveguide geometry. The pro- 
posed scheme relies mostly on the possibility of vanish- 
ing group velocity, i.e. diverging density of states, for 
the guided modes, at a frequency cutoff (or bandgap) of 
the waveguide. An important innovation of this work is 
the derivation of a non-perturbative analytic theory for 
RDDI near a cutoff of the photonic spectrum. The theory 
exhibits non-Markovian features, particulary population 
loss of the atoms by incomplete decay and the resulting 
reduction of entanglement, in agreement with numerical 
simulations. 

Possible manifestations of the predicted effect include 
high-concurrence entanglement as well as energy transfer 
between dipoles at giant separations. The analysis and 
the potential realizations discussed above suggest that 
the effect is significant for a wide range of atomic and 
waveguide parameters, constrained only by the tradeoff 



between interaction strength and the maximal achievable 
entanglement. 

We acknowledge the support of DIP, ISF and the Wolf- 
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APPENDIX 

A. Dipole-dipole interaction for arbitrary oriented 
dipoles 

In the main text we considered the case where the 
dipoles are oriented in the z direction. For a general 
orientation, we need to consider all the TE /TM mn ^ z 
modes with their normalized spatial functions [20] . 
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where A = ab is the transverse area of the waveguide. Inserting these mode functions into Eq. Q, we obtain the 
bath spectrum, 
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where Q(x) is the Heaviside step function. The effective dipole moments read 
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with dj = d ■ ej and x a , y a the transverse position of atom a. In order to find the RDDI A aa i 
we recall Eq. Q, and find by contour integration methods, 
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with £ mn = 
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B. Numerical simulations 

We performed direct numerical simulations of the 
Schrodinger equation for the Hamiltonian from Eq. (1), 
taking only the dominant TMn mode into account. The 
dipole couplings g k relate to the Id spectrum, from Eq. 
(7), by g UM — G aa {uj)duje lkzZa , where duj is the sam- 
pling resolution used to discretize the frequency space 
uj. The initial atomic state is \e\,g2) where the modes 
are in the vacuum |0). By taking the rotating wave ap- 
proximation |16j . i.e. neglecting non-energy-conserving 
Hamiltonian terms of the form a + d\ J1 a^d^, we restrict 
ourselves to the single-excitation Hilbert space, | e% , #2 > 0} , 



191:62,0) and {|<?i,<72j Vw}, which is solved numeri- 
cally. 



C. Validity of the Markov approximation 

The dissipative and dispersive coefficients, 7 QQ < and 
A aa >, can be obtained by deriving the master equation 
PT5l ITT] for the atoms' density matrix. Equivalently, here 
we will use instead the latter, second order perturbation 
theory for the transition amplitude. We begin with Eq. 
(23) on page 28 of Ref. [16], 
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where U^', is the second order contribution to the transition amplitude from the state where only atom a is excited 
to the state where only atom a' is excited, T is the interaction time, and 
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Recalling the definition of the bath spectrum in Eq. ([2]), we can rewrite W aa i as 



W aa > (w) = hm^_> + 



duj 



1 G aa '(oj') 



1 G aa i(u)) 



<jj — uj' — 17] J uj — 2ui a — uj 1 — it] 
Using the relation lim ?) _ > .Q+ ^jn- = iir5(x) + under integration, we obtain 



(19) 



(20) 



-i^7aa'M - i\p aa '{u - 2uj a ) - A aa > (u) - A aa ' {uj - 2uj a ) 



(21) 



with S t (uj) = J^ 2 /2 dte- tuJt being a sine function of width 1/T and amplitude T, and 

7ao'(^) = 2nG aa >(uj) ; A aa ,(uj) = P / dJ — ^ . 



(22) 
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In the limit T — > oo, i.e. <5t(w) ~ S(lj), we recover the Markovian results j aa > = 7 ac . (w ) and A aa > — A aa i(u] a ) + 
A aa >{— uj a ) [noting that G aa '(uj < 0) = 0]. Let us specify when such a limit is reasonable. Consider T as the time- 
resolution we are interested in, i.e. T is much smaller than the typical time-scale of the atomic dynamics. Nevertheless, 
we assume that T is sufficiently large, such that in a width 1/T of 8j.(oj — uj a ) around u> a , Jaa'(u), A aa >(u)) do not 
change appreciably. Then, we can expand "f aa i (w), A aa > (w) around ui a (and also around — ui a for A aa i) and get 

/ __■#(_> - u a )A aa , (_») oc A aa , (a;.) + 0[ aa T l a) ), (23) 



where a similar result is obtained for j aa r . For the 
Markovian approximation to be valid, we demand that 
the lowest order relative correction for the Markovian re- 
sult is small, 



KM 1 

Aaa'(Wa) T 2 



< 1. 



(24) 



As a typical atomic dynamics time-scale, for the case of 
RDDI, we may take A aa * . Then, using it in (24 1, we 
obtain the condition of validity in Eq. (12). 



D. Non-Markovian theory 

Taking the Laplace transform of the dynamical equa- 
tions for ai(t), a%{t) and i)fc(t) with the initial conditions 
oi (0) = 1, o 2 (0) = b k (0) = 0, we find 



oi(s) 



s + Jn(s) + iu a 



Jl2{s)J 2 l{s) 
S + J22(s) + lU a 



> (25) 



Here Oi(s) is the Laplace transform of ai(i) and 

Jaa'(s) — J2k 9k s +iu}k ' ^ e note that by virtue of 
Eq. Q, J aQ '(— ibj a ) = -iA aa ^_. As before, we con- 
sider the spectrum in Eq. ^ for m = l,n = 1. 
Since uj a is close to the cutoff un, the main contribu- 
tion to RDDI comes from frequencies near w u so that 
we approximat e the denominator of the spectrum by 
-^/(w/wii) 2 - la \/2yJw/u)ii — 1. After performing the 
integrals in J aQ '(s), using the approximated spectrum, 
we obtain 



ai(s) 
n{u) 
d{u) 



di(u) 



n{u) 
= dju) 

3 + W a u - 



2V2 



21F a u J 



V2 



V2 



(26) 



with it = y/^iy/a + ioJn, W a = u a — W11 and F(u) = 
( e -2(*i-;s 2 )(VwIT/c)u _ In order to p er f onn the inverse 

Laplace transform we first expand F(u) in a Taylor series: 
in order to still satisfy the Laplace initial value theorem, 



ax(t = + ) = lim.5_j.00 soti{s), the expansion is taken up 
to 5th order. Then, expanding 01(14) in partial fractions 



«i 



; Cj = c( Uj ) ; c(_) = (27) 



where Uy are the roots of d(u), and using the inverse 
transform of l/{\fs + a) [3D], we finally obtain 



01 (t) = v/ie"^ 11 ' ^ 9 



1 2 

— + v / iu J -e™3*erfc(— Viu 3 Vi) 
nt 

(28) 



E. Metal waveguide realization: ohmic losses 

We consider ohmic losses on the four conducting plates 
that make up the waveguide. The dissipated power per 
unit area of a plate is given by 



dP loss /dS = 0.5\Js\ 2 R s 



(29) 



where S is the area, R s its surface resistance [3T]. In 
order to find the surface current J s we should first find 
the electric field of the dipole inside the waveguide. As- 
suming, as before, that the dipole is oriented to the z 
direction, its field is a superposition of evanescent TM mn 
modes of a single cu a < uj mn photon, 



E?rin (*") 



HuJa 

2^ 



-u 



TM 



(30) 



where u^„ w (r) is given by Eq. (14 1 with k = 
( 1 /c) yj ui 2 m — uj 2 replacing — ik z and 2k replacing 1/L. 
We then find the corresponding magnetic field using the 
Maxwell equations for TM modes [20, 31J, 



-iw a e Vj_ x (E m „ • e z ), (31) 



H mn (r) 



Vi = d x e x + d y e y being the curl operator in the xy 
plane. The surface currents on the plates are found from 
the surface boundary conditions for the magnetic fields, 
J s = e n x H, with e„ the normal to the surface. Fi- 



nally we integrate Eq. ( |29[ ) over the plate area, e.g., 

2f™dztfdx0.5\J s \ 2 R s . 



for the plate at y = b, Pi OS s 



By defining the dissipation rate per TM mn mode as dispersion as 



mn 
I Loss 



mn 

I LOSS 



Pi oss /(fkj a ), we 


find for 


all four plates, 




2 




2 






L V n a J 


Moo 


i(£!) 2 + i. 


H a 



(32) 



Then, for the case a = b, the total dissipation of a single 
photon field from the atom is bounded by 2 . 



F. Fiber-Bragg-grating realization 

We briefly show how we can relate the fiber-Bragg- 
grating case to the theory derived for the MWG in the 
main text. The dispersion of a transverse fiber-mode with 
a Bragg-grating is [29]. 



ui{k z ) - uj b 



where ks = ^>b/c = ir/{kn) is the Bragg wavevector, A 
the grating period, n the average refractive index and An 
the index difference of the grating. Near the upper cutoff 
of the bandgap, k z is close to fee and we approximate the 




Lo(k z ) w lo u + B(k z - k B ] 2 



(34) 



where uj u = lob{^- + 0.5An/n) is the upper bandedge and 
B = (I ) 2 (sk) ZT- Then ' the density of states is 



duj 



4Ar 



V( w K) - 1 



(35) 



where ujb ~ w u was taken. There are three terms on 
the right-hand-side of Eq. ( 35 ) : the first is a linear dis- 



persion contribution of a mode with group velocity c/n, 
while the second increases the usual density of states by 
a constant factor. The third term is the divergence due 
to the bandedge. The spectrum of the fiber mode will 
then have the form [see Eq. (7)] 



G aa (uj) 



1 



2tt v^MO-r 



(36) 



where T u is similar to the free space spontaneous emission 
rate. This is the spectrum assumed in our non-Markovian 
theory for the MWG, with uj u ,T u replacing wu,rn. 
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